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ABSTRACT 

We present a long (116 ks) Suzaku observation of the Broad-Line Radio 
Galaxy (BLRG) 3C 382 acquired in April 2007. A Swift BAT spectrum in 15- 
200 keV from the 58-month survey is also analyzed, together with an archival 
XMM-Newton EPIC exposure of 20 ks obtained one year after Suzaku. Our main 
result is the finding with Suzaku of a broad Fe K line with a relativistic pro- 
file consistent with emission from an accretion disk at tens of gravitational radii 
from the central black hole. The XIS data indicate emission from highly ionized 
iron and allow us to set tight, albeit model-dependent, constraints on the inner 
and outer radii of the disk reflecting region, rj„ ~ lOr^ and Vout — 20rg, respec- 
tively, and on the disk inclination, i ~ 30°. Two ionized reflection components 
are possibly observed, with similar contributions of ~10% to the total contin- 
uum. A highly ionized one, with log^~3 erg s~^ cm, which successfully models 
the relativistic line and a mildly ionized one, with log^~1.5 erg s~^ cm, which 
models the narrow Fe Ka and high energy hump. When both these components 
are included, there is no further requirement for an additional black body soft 
excess below 2 keV. The Suzaku data confirm the presence of a warm absorber 
previously known from grating studies. After accounting for all the spectral fea- 
tures, the intrinsic photon index of the X-ray continuum is Fa, ~ 1.8 with a cutoff 
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energy at ~ 200 keV, consistent with Comptonization models and excluding jet- 
related emission up to these energies. Comparison of the X-ray properties of 
3C 382 and other BLRGs to Seyferts recently observed with Suzaku and BAT 
confirms the idea that the distinction between radio-loud and radio-quiet AGN 
at X-rays is blurred. The two classes form a continuum distribution in terms 
of X-ray photon index, reflection strength, and Fe K line width (related to the 
disk emission radius), with BLRGs clustered at one end of the distribution. This 
points to a common structure of the central engine, with only a few fundamental 
parameter (s) responsible for the radio- loud/radio-quiet division. The black hole 
spin, and in particular its rotation compared to the disk's, may be a key one. 

Subject Headings: Galaxies: active — galaxies: radio - galaxies: individual — X-rays: 
galaxies 

1. Introduction 

It is by now widely accepted that Active Galactic Nuclei (AGN), including radio-loud 
sources, are powered by accretion of the host galaxy gas onto a central super-massive black 
hole (e.g., Blandford 1985). The rich variety of AGN phenomenology has been explained 
in terms of orientation with respect to the axis of an obscuring equatorial thick matter 
surrounding the nucleus (Antonucci 1993; Urry & Padovani 1995), yielding type-1 sources 
for face-on views and type-2 objects for edge-on views. In addition, in radio-loud sources a 
relativistic jet is present, connecting the innermost regions near the black hole to the radio 
galaxy perifery. The jet angle - defined as the angle between the jet axis and the line of sight 
- increases from blazars to Broad- and Narrow-Line Radio Galaxies, roughly corresponding 
to an increase of obscuration degree. As the jet angle increases, the importance of its emission 
decreases, due to beaming effects. 

X-ray spectroscopy is an effective tool to investigate the central engines of AGN. In- 
deed, previous X-ray observations of BLRGs in the 1990s with ASCA, RXTE, and BeppoSAX 
established that these sources exhibit Seyfert-like spectra with subtle but significant differ- 
ences, i.e., flatter X-ray continua and weaker reflection features than in radio-quiet (see 
Sambruna, Eracleous, & Mushotzky 2002; Ballantyne 2007, and references therein, for a re- 
view). In particular, the Fe Ka emission line around ~6.4 keV was observed to be narrower 
and of lower EW than in Seyferts, and generally unresolved at the limited resolution of the 
pie- Suzaku detectors. No clear evidence for a relativistic accretion disk profile was found for 
the Fe Ka line in any BLRGs, or other RL AGN. These results raised the currently held 
scenario that some intrinsic, fundamental difference must exist in the structure of the accre- 
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tion flow between the two classes of AGN, with RQ being dominated by standard optically 
thick, geometrically thin disks and RL having radiatively inefficient, ADAF-like, inner disks 
(see discussion in Eracleous et al. 2000; Ballantyne 2007). Alternatively or in addition, a 
non-thermal jet contribution was also possible. 

The advent of Suzaku with its broad-band coverage in 0.3-100 keV and improved sen- 
sitivity especially in the critical Fe K region, 6-7 keV, has provided a golden opportunity 
to deepen our understanding of the central engines of RL AGN. Our group has secured GO 
observations of all the bright, nearby {z < 0.1) BLRGs with relatively deep, 100 ks, Suzaku 
exposures. The Suzaku observations of our program targets arc described in separate pa- 
pers that concentrate on one source at the time, as each classical BLRG presents somewhat 
unique X-ray and multi-wavelength properties (3c 390.3, Sambruna et al. 2009, S09 here- 
after; 3C 445, Braito et al. 2011; Reeves et al. 2010; 3C 111, Ballo et al. in prep.). In 
this paper, we focus on 3C 382, a giant lobe-dominated radio galaxy well-known for its soft 
excess and rapid X-ray variability from previous observations (sec § 2). Besides its Suzaku 
data, we also discuss its XMM-Newton EPIC data and a Swift BAT spectrum (see below), 
both unpublished. 

The Suzaku observation of 3C 382 confirms the complexity of its spectrum and its simi- 
larity to Seyferts. Remarkably, we find compelling evidence in the XIS data for a relativistic 
Fe Ka line profile, with excellent constraints on the disk inclination and inner/outer radii 
(§ 5.3). This is the first time that a such a detailed profile is detected in a BLRG, and indeed 
in a radio-loud AGN. The implications for the structure of the central engine and models of 
jet formation are far reaching (§8). 

The paper is organized as follows. After describing the source properties and previous 
observations in § 2, in § 3 the data reduction of the new observations is presented. The Suzaku 
timing analysis is reported in § 4. We describe the spectral fits to the various datasets in 
§ 5 for the combined Suzaku and Swift BAT and in § 6 for XMM-Newton. In § 7 we report 
the summary of the results and the discussion following in § 8. Finally, in Appendix A we 
report the background and cahbration tests for the XIS cameras. Throughout this paper, a 
concordance cosmology with Hq = 71 km s"^ Mpc""^, Qa=0.73, and fij„=0.27 (Spergel et al. 
2003) is adopted. The energy spectral index, a, is defined such that oc The photon 
index is F = a -|- 1. 
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2. 3C 382 and Previous X-ray Observations 

The lobe-dominated, Fanaroff-Riley II radio galaxy 3C 382 {z = 0.0579) has a double- 
lobe structure, with a clear jet in the northern lobe ending in a hotspot. While a counter-jet is 
not clear, a hotspot is also detected in the southern lobe, with a total size between hotspots at 
3.89 GHz of 179" (Hardcastle et al. 1998). Prom this size, a lower limit to the jet inclination 
of ^ = 15° is inferred (Eracleous & Halpern 1998). While it is common to assume that in 
BLRGs the jet axis is parallel to the disk axis, so that 9 ^ i, with i the disk inclination with 
respect to the line of sight, in this paper we will leave i free to vary during the spectral fits 
to the X-ray data whenever possible. 

Optically, 3C 382 is identified with a disturbed elliptical galaxy dominated by a very 
bright and unresolved nucleus (Matthews, Morgan Sz Schmidt 1964; Martel et al. 1999), 
located in a moderately rich environment (Longair & Seldner 1979). The optical spectra 
show a strong continuum and prominent broad lines photo-ionized by a power-law type of 
spectrum (Saunders et al. 1989; Tadhunter, Fosbury & Quinn 1989), with FWHM ~ 12,000 
km/s for the Ha line (Eracleous & Halpern 1994). A recent estimate using the luminosity 
of the host galaxy places the mass of the central black hole around Mbh ~ 10^ M© within 
40% (Marchesini et al. 2004). 

3C 382 is a variable source at X-ray (Dower et al. 1980; Barr & Giommi 1992), radio 
(Strom, Willis & Willis 1978), optical, and UV frequencies (Puschell 1981; Tadhunter, Perez 
& Fosbury 1986). At X-rays it is well studied, and was observed by all the major X-ray 
observatories before Suzaku. Flux and spectral variability is observed on time-scales shorter 
than a day, with a trend of spectral softening for increasing intensity (Gliozzi et al. 2007). 

Previous low-sensitivity observations of 3C 382 indicated that its X-ray spectrum was 
remarkably similar to Seyferts (Reynolds 1997). Aside from the variable soft excess below 
2 keV, 3C 382 exhibits a warm absorber with ionization parameter log^ ~ 2.5 erg cm s~^, 
column density ~ 1 — 3 x 10^^ cm~^, and outfiow velocity ~1,000 km/s, clearly detected 
in recent high-resolution dispersion spectra (Reeves et al. 2009; Torresi et al. 2010). The 
X-ray continuum can be described by a power law with the "canonical" photon index P 1.8 
(e.g., Sambruna et al. 1999, S99 hereafter) and weak reflection features, including a Compton 
hump above 10 keV with reflection fraction R ~ 0.2 — 0.6 (Gliozzi et al. 2007) and an Fe 
Ka line. The profile of the Fe K line has remained, so far, ambiguous due to its inherent 
complexity, as we show here (§ 5.1), and the poor sensitivity of the pre-Suzaku detectors. 
For example, using ASCA an unusually broad (Gaussian width (Tg~1 keV) line was inferred 
(S99; Reynolds 1997) while only a narrow component with (Tg~100 eV was detected with 
RXTE, BeppoSAX, and the Chandra HETGS (Eracleous et al. 2000; Grandi et al. 2001; 
Gfiozzi et al. 2007). 



-5 - 



The source exhibits a strong, variable X-ray excess over the extrapolation of the power- 
law continuum from the hard X-rays below 2 keV. Extended soft X-ray emission was detected 
with the ROSATRRl in 0.3-2.4 keV (Prieto 2000), later confirmed with Chandra (Gliozzi et 
al. 2007). The contribution to the spectral soft excess from the extended emission, however, 
is negligible as shown by the point-like ACIS-S image in 0.5-2 keV (Gliozzi et al. 2007). 

Finally, while other BLRGs have been detected at GeV gamma-rays (3C 111; Hartman 
et al. 2008), no signal has been reported so far from 3C 382 with Fermi (Abdo et al. 2010). 
The lack of gamma-ray flux and the double-lobe radio structure are strong indicators that 
the jet in this source is misaligned and most likely does not contaminate the nuclear emission, 
as supported by the Seyfert-like X-ray spectrum from the ASCA and RXTE era. 

3. New X-ray Observations: Data Reduction 

The log of the X-ray observations is reported in Table 1. The exposures are after 
data screening and the source count rates after background subtraction, according to the 
procedures described below. 

3.1. Suzaku XIS 

Suzaku observed 3C 382 on April 27th, 2007 for a net exposure time after screening of 
~116 ks. We used the cleaned event files obtained from version 2 of the Suzaku pipeline 
processing. Standard screening criteria were used, namely, only events outside the South 
Atlantic Anomaly (SAA) as well as with an Earth elevation angle (ELV) > 5° were retained, 
and Earth day-time elevation angles (DYE_ELV) > 20°. Furthermore, data within 256 s of 
the SAA passage were excluded and a cut-off rigidity of > 6 GV was adopted. The main 
parameters of the XIS observations are reported in Table 1. 

The XIS spectra of 3C 382 were extracted from a circular region of 2.9' radius centered 
on the source. Background spectra were extracted from four circular regions offset from the 
main target and avoiding the serendipitous source and the calibration sources. The combined 
area of these four background regions is twice the area of the main target region. The XIS 
response (rmf) and ancillary response (arf) files were produced, using the latest cahbration 
files available, with the ftools tasks xisrmfgen and xissimarfgen, respectively The source 
spectra from the front illuminated (FI) CCDs were summed, and fitted jointly with the back 
illuminated (Bl), the XISl, spectrum after verifying that the data from the FI and BI CCDs 
were consistent with each other. In fact, we checked that the 2-10 keV continuum slope and 
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fluxes in the two cases are consistent within the 2%. The net XIS background-subtracted 
source spectra were binned with a minimum of 25 counts per energy bin to ensure that the 
goodness of fit can be applied to the spectral analysis. This binning is maintained throughout 
the spectral analysis and only in some figures it was increased for clarity, when specified. 
Consistent results were found by varying the binning in a range up to 100 counts/bin. 

3.2. Suzaku HXD 

For the HXD-PIN data reduction and analysis we followed the latest Suzaku data reduc- 
tion guide (the ABC guide Version 2)0. We used the rev2 data, which include all four cluster 
units, and the best background available, which account for the instrumental background 
(Kokubun et al. 2007)@. The source and background spectra were extracted within the com- 
mon good time interval and the source spectrum was corrected for the detector dead-time. 
The net exposure time after screening was 113 ks. 

The contribution of the diffuse cosmic X-ray background counts was simulated using the 
spectral form of Boldt (1987), assuming the response matrix for diffuse emission, and then 
added to the instrumental one. Two instrumental background files are available; background 
A or "quick" background, and background D or "tuned" background. We adopted the latter 
which is the latest release and which suffers lower systematic uncertainties of about 1.3%, 
corresponding to about half the uncertainty of the first release of the Non X-ray Background. 
With this choice of background, 3C 382 is detected up to 70 keV with the PIN at a level of 
~20% above the background. The count rate in 10-30 keV is 0.135 ± 0.003 c/s. The HXD 
PIN spectrum was binned in order to have a signal-to-noise ratio greater then 10 in each 
bin, and the latest response file released by the instrumental team was used. 

For a comparison with the Swift BAT spectrum, we considered the PIN spectrum in 
the overlapping energy band of the two detectors, E=15-70 keV (see § 3.4). In this band, 
the PIN spectrum can be fitted with a single power law with Fp/Ar = 1-95 ± 0.14, yielding a 
flux of 6.37l5'23 ^ 10~^^ erg cm"^ s~^, with 90% errors. 



^http://heasarc. gsfc.nasa.gov/docs/suzaku/analysis/abc/ 

^ftp:/ /legacy.gsfc.nasa.gov/suzaku/doc/hxd/suzakumemo-2008-03.pdf 
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3.3. XMM-Newton 

XMM-Newton observed 3C 382 on April 28th, 2008, for a net exposure after screening of 
20 ks. The analysis of the RGS spectrum has already been performed by Torresi et al. (2010) 
and is not reported here. 

We downloaded and reduced the data of the three EPIC cameras, pn, MOSl and M0S2. 
The data reduction was performed following the standard procedure with the XMM-SAS v. 
8.0.1 package. We checked the observation for high background contamination, looking for 
flares in the light curves at energies greater than 10 keV. We excluded these bad intervals 
from the successive analysis, and extracted the source photons from a circular region of 40", 
while the background ones were collected from an adjacent source free circular region of the 
same size. Only single and double events were selected. 

Using the SAS task epatplot we checked that the pile-up fraction of the EPIC pn is 
negligible, less than 1%. Instead, given the X-ray brightness of the source, the pile-up 
fraction of the EPIC MOS observations is significant, ~50%. The main result of pile-up 
is to introduce an overall distortion of the observed continuum, as it basically consists in 
counting two or more low-energy photons as a single high-energy photon. In fact, we find 
that the 2-10 keV power-law continuum of the pn has a r~1.74, while for the MOS r~1.6. 
Therefore, given the significant pileup of the MOS and and the ^25% lower S/N ratio of 
the MOS data, due to the short exposure and lower effective area, we limited the analysis of 
XMM-Newton only to data from the EPIC pn camera. 

Similarly to the XIS, the net EPIC pn background-subtracted source spectrum was 
binned with a minimum of 25 counts per energy bin to ensure that the goodness of fit 
can be applied to the spectral analysis This binning is maintained throughout the successive 
spectral analysis and increased only in some figures for clarity of presentation, when specified. 
Consistent results were found by varying the binning in a range up to 100 counts/bin. 



3.4. Swift BAT Observations 

The BAT spectrum was derived from the 58-month hard X-ray surve}0. The data 
reduction and extraction procedure of the 8-channel spectrum is described in Baumgartner et 
al. (2010). To fit the BAT spectrum, we used the latest calibration response diagonal_8 . rsp 
and background files as of December 2010. 



^http://heasarc. gsfc.nasa.gov/docs/swift/results/bs58nion/ 
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The source is clearly detected in the energy band E=15-200 keV with a S/N=35 and 
a count rate (1.35 ± 0.05) x 10~^ c/s. A fit to the mean 58-month BAT spectrum in the 
15-70 keV band with a single power law model yields a photon index Tbat = 2. 12 ±0.20 and 
a flux of 5.10~^^'ll X 10^^-*^ erg cm^^ s^-*^, at 90% uncertainty. These values are consistent with 
those derived for the Suzaku HXD PIN in the same energy range (see §3.2), justifying joint 
fits using both instruments. This is especially important in view of the extended sensitivity 
to higher energies provided by the BAT, which will allow us to set better constraints on 
reflection models. 



4. Timing Analysis 

Figure 1 shows the Suzaku XlSO-3 light curves of 3C 382 at soft and hard X-rays, 
respectively, with the flux in bins of 5760 sec (1 orbit). The data from the three detectors 
are plotted on separate panels for comparison. Weak flux variations are seen in both energy 
bands for all three detectors, with small- amplitude variations of the order of few percent. 
When the ratio of the counts in the two energy bands is plotted against time or total count 
rate, there is no evidence for spectral variability. We note the presence of two possible short 
time-scale flares of ~5% amplitude in the XIS 3. We ascribe them to spurious instrumental 
fluctuations. However, as discussed in the Appendix, they have a negligible influence on 
the spectral analysis of the average spectrum reported here. Aside from the small apparent 
flares, the overall shape of the light curves of the different XIS detectors is consistent, with 
each showing a gradual increase in counts over the Suzaku observation. 

The variability analysis was repeated for the XMM-Newton EPIC data, despite the 
short exposure affected by the large background ffares. No evidence for ffux changes was 
found. Table 2 reports the soft (0.5-2 keV) and hard (2-10 keV) fluxes of 3C 382 from 
the Suzaku and XMM-Newton observations. The reported fluxes are both observed and 
corrected for any intervening absorption and were derived from the best-fit models to the 
0.6-200 keV Suzaku+BAT and 0.5-10 keV EPIC pn spectra (see § 5.4 and § 6). Observed 
and absorption corrected luminosities are reported in Table 2 as well. Comparing both fiuxes 
and luminosities in all energy bands, it is clear that the source was in similar intensity states 
at the epoch of the Suzaku and XMM-Newton observations. 
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5. The Suzaku Broad-Band Spectrum 

The lack of significant spectral variability derived from a model-independent timing 
analysis of the XIS data (§ 4) indicates that the spectral properties of 3C 382 can be safely 
investigated using a time-averaged spectrum accumulated over the entire 116 ks exposure. 
This spectrum provides the largest signal-to-noise ratio available from this observation, al- 
lowing the most detailed spectral analysis of the Fe K line region and broad-band continuum, 
the main goals of our project. 

All spectral fits were performed using XSPEC v.12.6.0. The significance of adding free 
parameters to the model was evaluated with the F-test. We included only model components 
with a statistical probability Pf>99%. All uncertainties quoted are la for one parameter 
of interest, unless otherwise stated. The energy of the lines is reported in the source rest- 
frame, if not otherwise stated. The data from the XIS CCDs were fitted jointly, excluding 
the energy range E=1.5-2 keV, around the Si K edge, which is known to be affected by 
calibration issues (Koyama et al. 2007). When fitting the broad-band spectrum, a constant 
factor was introduced in the model to take into account the cross-normalization between the 
different instruments. In particular, a constant factor of ~ 1.16 was adopted for the XIS/PIN 
and of ~0.9 for the PIN/BAT. Throughout the spectral analysis we always include a Galactic 
absorption towards 3C 382 of Nh^7A x 10^° cm'^ (e.g.. Dickey & Lockman 1990; Kalberla 
et al. 2005). 

The plan of the analysis is as follows. First, we consider the combined broad-band 
Suzaku and Swift BAT spectra in the E=0. 6-200 keV, excluding the Fe K band, and derive 
an estimate of the underlying X-ray continuum. Then, we consider also the Fe K band and 
perform a search for additional emission/absorption features in the spectrum. Subsequently, 
we perform fits to the full broad-band spectrum with a physical model including all the 
spectral components identified in the previous steps plus an additional relativistic line and 
smeared ionized refiection from the putative accretion disk. Finally, we test the simultaneous 
alternative modeling of the soft excess, Fe K lines and cold refiection with a second lowly 
ionized refiection component. 

5.1. The Continuum Determination 

The combination of the Suzaku and Swift BAT instruments provides an unprecedentedly 
wide and sensitive spectral coverage from 0.6 keV up to 200 keV. This is fundamental for a 
correct determination of the underlying X-ray continuum. From this analysis, we exclude the 
Fe K energy range, E=5.5-7.5 keV, because, as discussed in the next section, the presence 
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of spectral complexities in this interval can affect the best-fit continuum model. 

In the upper panel of Figure 2 we show the ratio of the combined XIS03, XISl, PIN, 
and BAT spectra with respect to a Galactic absorbed power-law continuum. As seen, this 
simple model does not provide a sufficient parameterization of the data (x^/do/=5027/3166) 
and strong deviations are present at both low and high energies. The resultant power-law 
photon index is r~1.88. 

At energies lower than ~2 keV there is evidence for a soft excess component. This 
is expected also from previous studies of 3C 382, see § 2. Following previous authors, we 
initially model the soft excess with a phenomenological blackbody with best-fit temperature 
A;T~100 eV. This is highly required, with a Ax^~1407 for two additional model parameters. 
The alternative phenomenological modeling with a power-law provides a comparable fit. 
However, the physical interpretation of the resultant very steep photon index of r~5 is more 
ambiguous and therefore in the following we keep the blackbody parameterization. 

The excess flux in the upper panel of Figure 2 between 10-40 keV is instead indicative of 
a reflection component. This is also required on a physical basis by the presence of a narrow 
Fe Ka line (e.g., Lightman et al. 1988), as discussed in the next section. Therefore, we added 
a cold reflection component using pexrav in XSPEC. We assume an inclination angle of i=30°, 
typical for BLRGs, and standard solar abundances. This component is highly required, with 
a Ax^~66 for two model parameters, significant at >99.99%. The power-law continuum 
slope is now r~1.85, the high energy cut-off Ec~200-300 keV, and the refiection fraction 
from pexrav is R~0.36. The final fit is good, with a reduced chi-squared very close to unity, 
X^/(io/=3554/3162. The ratios of the spectra with respect to this best-fit continuum model, 
including the Fe K band, are shown in the middle panel of Figure 2. 

In conclusion, the broad-band X-ray continuum of 3C 382 is well parameterized by 
a Galactic absorbed power-law, plus a blackbody soft excess and a reflection component. 
We use this as the baseline continuum model for the successive search of additional emis- 
sion/absorption spectral features. 

5.2. Search for Emission/ Absorption Spectral Features 

Looking at the spectral ratios in the middle panel of Figure 2, weak dips indicating 

absorption can be seen at E<2 keV, as well as further spectral complexities in the Fe K 
region. The BLRG 3C 382 is known to exhibit absorption features at energies below 2 keV 
due to a warm absorber (Reynolds 1997), recently studied with the Chandra and XMM- 
Newton gratings (Reeves et al. 2009; Torresi et al. 2010). Despite the XIS lower resolution 
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compared to the gratings, absorption features are visible in Figure 2 (middle panel) around 
0.8-1 keV, possibly due to ionized O and Fe L. We thus include a warm absorber component 
in our fit, modeled by an XSTAR grid (Kallman & Bautista 2001). Given the limited S/N 
and energy resolution available with respect to the gratings studies, we fix the gas outflow 
velocity to 1,000 km/s, consistent with what found by Reeves et al. (2009) and Torresi et 
al. (2010). We assumed an ionizing X-ray power- law with photon index r=2 and turbulence 
velocity of 500 km/s, similar to what assumed by Reeves et al. (2009). 

The inclusion of the warm absorber component is highly significant, Ax^=30 for two 
additional parameters, corresponding to a detection probability >99.99%. Wc estimate a 
column density N//~6 x 10^° cm^'^ and an ionization parameter log.^~2.44 erg s^^ cm. 

We now reintroduce the E=5. 5-7.5 keV energy band in the fits. Up to now the fit 
statistics is x^/dof = 3524/3160. Including the Fc K band, this changes to x^/dof = 
4760/4177. The middle panel of Figure 2 and upper panel of Figure 3 show the presence of 
spectral complexities in the Fe K band. In particular, a prominent emission line is visible at 
rest-frame energy ~6.4 keV, together with several other emission lines in the energy range 
~5.5-7.5 keV. 

As a first step, we modeled these emission features with Gaussians. The most promi- 
nent line at E= 6.41 ± 0.01 keV is clearly identified with the Fe Ka from lowly ionized, 
fiuorescent material. The line is resolved, with a width of (7(3=1 18I13 ^V, equivalent width 
EW=68 ± 9 eV, and fiux I=(3.3 ± 0.2) x 10"^ ph s"^ cm-^. The inclusion of this feature 
yields a Ax^=233 for three additional model parameters, which corresponds to a detection 
probability >99.99%. 

Next, we added the associated Fe K/3 component. In the fit, the energy was fixed to the 
expected value of E=7.06 keV, the width was set equal to that of the Fe Ka and its intensity 
to ~13% of the Fe Ka (e.g., Molendi, Bianchi & Matt 2003 and references therein). The 
addition of the Fe K/3 improves the fit by Ax^=13. 

Three additional emission features are still present in the spectral residuals at energies of 
E~5.8, 6.9, and 7.5 keV. To account for these, we added three more Gaussians to the model. 

The two lines at higher energies are not resolved. Their energies arc at E=6.91 ± 0.02 keV 
and E=7.51 ± 0.03 kcV and have equivalent widths of 17 ± 6 eV and 18 ± 7 cV. respectively. 
The Ax^ is 21 and 19, for two additional model parameters, which correspond to detection 
probabilities of ^99.99%. The lower energy line at E=5.70±0.09 keV is instead broader and 
resolved, with a width of a= 385^57 eV and equivalent width of 55 ± 12 eV. Its significance 
is higher, Ax^=58 for three additional parameters, corresponding to a detection probability 
>99.99%. 
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The fit with the five Gaussian fines is purely phenomenological, in order to have a first 
representation of the Fe K complex and a determination of the lines parameters and their 
significance. All the Gaussian lines are detected with a probability >99.99%. Goncerning 
the identification of the three additional emission lines, the feature at 7.5 keV may well 
be identified with emission from lowly ionized Ni Ka, whose atomic transition energy is 
expected at E~7.5 keV (e.g., Molendi, Bianchi & Matt 2003 and references therein). The 
two remaining fines at 6.9 keV and 5.8 keV are at energies not directly consistent with any 
significant atomic transitions and therefore have a less obvious physical interpretation. The 
closest atomic fines to the 6.9 keV line are the Fe XXV Hea (ls^-ls2p) at E~6.7 keV and 
Fe XXVI Lya (ls-2p) at E~6.97 keV (KaUman et al. 2004). However, the identification 
with one of these lines would require a blue/red-shift of ~0.03c or ~0.01c, respectively. The 
5.8 keV line is even more mysterious, being out of the range expected for Fe K lines, from 
6.4 keV up to 7 keV, depending on the ionization. If identified with Fe K emission, the large 
width and red-shift could suggest strong relativistic effects and an origin close to the central 
black hole (see next section). 

The inclusion of the warm absorber and the emission lines in the Fe K band on the 
continuum model provides a very good fit, with a x^/(io/=4416/4167. Therefore, the final 
best-fit model of the combined broad-band Suzaku and Swift BAT spectrum of 3G 382 is 
rather complex and is composed by a Galactic absorber power-law continuum, a black body 
soft excess, a cold refiection component, a warm absorber and five Gaussian emission lines 
in the Fe K band. The best-fit parameters are reported in Table 3. The spectral ratios with 
respect to this model are shown in the lower panel of Fiure. 2. 

5.3. A Relativistic Fe K Emission Line in 3C 382? 

Inspection of Figure 3 (upper panel) suggests the intriguing possibility that the 5.7 keV 
and 6.9 keV emission features represent two components of a single, broader feature, i.e. 
a relativistic Fe K line. More precisely, they would correspond to the broad red-wing and 
sharp blue-peak components of such a line, respectively (e.g., Fabian et al. 1989). As an 
initial test to this hypothesis wc replaced these two features with a single relativistic line 
component, diskline in XSPEC (Fabian et al. 1989), in the model discussed in the previous 
section and reported in Table 3. Therefore, now the Fe K band is composed by three narrow 
emission lines, an Fc Ka and associated Fe K/3 and Ni Ka, and a broad relativistic line (see 
middle panel of Figure 3). 

The diskline model considers a non-rotating Schwarzschild black hole, with spin param- 
eter a=0. Given that the fit is not sensitive to the power-law emissivity profile, we assumed 
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the typical value of 2.5 (e.g., Nandra et al. 2007; de la Calle Perez et al. 2010; Patrick 
et al. 2010). However, we checked that a different choice of f3 in the wide range from —5 
to does not significanly change the parameters estimates reported below, being always 
consistent at the 90% level. 

Interestingly, if the energy of the line is left free to vary, the fit yields a value larger 
than the expected 6.4 kcV for neutral Fe Ka emission. The formal best-fit ccntroid is at 
E=6.69 keV, ranging between 6.55 and 7.10 keV at 90% confidence. This indicates that the 
relativistic Fc K emission originates from ionized material, preferentially Fe XXV Hea at 
6.7 kcV rcst-framc, but wc can not clearly exclude also a contribution from Fe XXVI Lya at 
6.97 kcV. This supports earlier models claiming that the accretion disk in BLRGs is ionized 
(Ballantyne et al. 2002), and adds to previous evidence from the Suzaku observations of two 
other classical BLRGs, 3C 390.3 (Sambruna et al. 2009) and 3C 120 (Kataoka et al. 2007). 

Thus, we performed fits assuming, in turn, either emission from Fe XXV or Fe XXVI. In 
both cases the fit is highly improved, Ax^=60 for 4 additional parameters, corresponding to 
a detection significance >99.99%. The middle panel of Figure 3 illustrates the best-fit model 
for the three Gaussians and the diskline profile. Assuming Fe XXV (Fe XXVI) emission, we 
obtain good constraints on the disk inchnation i — 30° ± 1° (26° ± 2°), on the inner and 
outer radii of the emitting region. Tin = 12 ± 2 (10 ± 2 Vg) and rout = 23 ± 3 (18 ±3 Vg). 
Indeed, a tight estimate of i was expected from the 6.9 keV peak sharp fall to the blue (see 
Figure 3, middle panel), while the limited extension of the red "wing" indicates that the 
inner radius can not extend too further in close to the black hole. The equivalent width of 
the line is EW=80 ± 25 eV and the intensity I=(3.2 ± 0.4) x 10"^ ph s"^ cm"^^ 

We stress that the inclusion of the warm absorber has a negligible effect on the Fc K 
parameters (e.g., Reeves et al. 2004 and references therein). This can be explained by the 
low level of ionization and small column density of the warm absorber, which implies that 
only the light elements (C, O, Mg, Ne) are contributing to the absorption at energies below 
2 keV. 

The final fit with the relativistic line profile is very good, x^/dof =4435/4168 and the 
data/model ratio with respect to this model are shown in the lower panel of Figure 3. Given 
the high quality of the XIS data in 4-10 keV and good determination of the underlying 
continuum in the wide 0.6-200 keV band, we regard the detection of a broad, relativistic 
Fe K line in 3C 382 as robust. In support of this, we performed extensive calibration and 
background tests, which are described in the Appendix. 
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5.4. Detailed Modeling of the Relativistic Line and Ionized Reflection 

Here we discuss in more details the modeling of the broad Fe K line and possible ionized 
reflection component in the broad-band 0.6-200 keV spectrum. We first examine the broad 
Fe Ka profile in light of the possibility that it could be due to a rotating (either prograde or 
retrograde) black hole. This is possible by using the relativistic line code relline of Dauser 
et al. (2010), which allows also the black hole spin a to vary between —0.998 and +0.998. 
We use the final best-fit model discussed in § 5.2 and reported in Table 3, substituting the 
5.7 keV and 6.9 keV emission hnes with a single broad relativistic hne. 

The main parameters of relline are the rest-frame energy of the line, E, the power-law 
emissivity profile, (3, the black hole spin, a, the inner and outer radii of the reflecting region, 
rin and rout, the disk inclination, i, and the intensity of the emission line, /. At first, all 
the above parameters were free to vary during the fit. Only rj„, Tout, h I were constrained. 
We explored the parameter space in more detail, checking the effects of assuming different 
values. However, given the limited S/N of the present data and the complexity of the model, 
we are unable to constrain the spin, emissivity and energy all together. Nevertheless, good 
fits are obtained for |/3|~2-3. The line energy is E~6.6-7 keV at 90%, confirming ionized 
emission. Similar fits are obtained assuming Fe XXV Hea at E=6.7 keV or Fe XXVI Lya 
at E=6.97 keV, with shghtly smaller r^^ and ro„t in the latter case. 

We also explored the dependence of the line shape on the black hole spin by imposing 
that the inner radius coincide with the Innermost Stable Circular Orbit (ISCO) of the disk. 
This corresponds to assuming that the beginning of the reflection region is coincident with 
the inner radius of the disk and that the latter extends down to the ISCO. The ISCO 
assumes the following values depending on a: rjsco—^-'^^fg: Sr^ and 9rg for a=+0.998, 
and —0.998, respectively. A series of fits were performed with the inner radius linked to the 
ISCO parameter in the relline model in XSPEC, and leaving a free to vary. However, also 
in this case we can not still constrain all the parameters simultaneously and find a clear deep 
minimum in the statistical distribution. For instance, given the values in the intervals /3=1- 
3, E=6.7-6.97 keV and a=— 0.998 — 1-0.998, there are several minima, all within a difference 
of a few in A^^, and the fit does not reach a stable convergence. 

Therefore, despite the good quality of the present Suzaku data and the broad-band 
coverage up to 200 keV, no sensitive constraint was possible on the value of the black hole 
spin, i.e., we obtain similar fits for a = —0.998, a = 0, and a = +0.998. Detailed simulations 
show that much larger exposures are needed, ~ 300 ks, to achieve the goal of measuring a 
with sufficient accuracy to distinguish among a static, prograde or retrograde spinning hole. 

Given the results described so far, in the remaining fits to the 0.6-200 keV spectrum 
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of 3C 382 the broad Fe K line energy was fixed to E=6.7 keV, a—0, and ^ — —2.5, an 
intermediate value consistent with Seyferts (Nandra et al. 1997, Patrick ct al. 2010, de 
la Calle Perez et al. 2010). Including this component, the fit is improved by Ax^=60 for 
4 additional parameters, corresponding to a high detection probability >99.99%. This fit 
yields ri„=12 ±2 Vg and rout=23tt Tg, with t=30° ± 1°. The line intensity is J=(3.0 ± 0.4) x 
10^^ ph s^^ cm~^ and EW=76 ± 18 eV. This indicates that the broad, ionized Fe K line is 
produced from a narrow annulus of rin^lOvg and rout~20rg, i.e., at some distance from the 
black hole. We checked that a different choice of the /3 in the range —2 to —3 does not change 
significantly the parameters estimates, being always consistent within the 90% errors. 

As discussed above, the presence of ionized Fe K emission lines in the energy range 
6-7 keV provides strong motivation for including an additional ionized reflection component 
to our baseline continuum model. Sophisticated models are known in the literature that 
include the proper number of emission lines expected for the fitted ionization parameter, 
^ (e.g., Ross, Fabian & Young 1999; Ross & Fabian 2005; Garcia & Kallman 2010). Two 
such models, reflionx (Ross & Fabian 2005) and the more recent xillver (Garcia & Kallman 
2010), are available in XSPEC for detailed fits. We choose to use the latter, convolved with 
the relconv model of Dauser et al. (2010), which takes into account line blurring due to 
relativistic effects. We use a xillver table with log(^=0.8-3.8 erg s~^ cm, standard solar 
abundances, and an incident ionizing power-law with r=2. Gonsistent results were obtained 
using reflionx. 

The free parameters of the relconv component arc the power-law cmissivity profile, 
inner/outer radii on the disk surface, spin, and inclination angle. The free parameters 
of xillver are instead the ionization level of the material and the normalization. During 
the first fits to the 0.6-200 keV spectrum, these parameters were left free to vary. The fits 
yielded ri„~10 r^, rout— 25 T^g, inclination i~30°, and ionization parameter log^~3 erg s~^ cm, 
while the disk emissivity (3 was less constrained and the spin a was unconstrained. After 
performing a series of tests, we fixed /3=— 2.5, as no significant changes in the line and 
continuum parameters were observed within 90% uncertainties if the emissivity profile was 
left free to vary in the most probable range |/3|=2-3. We then linked rin=rjsc'0 in relconv, 
and performed three separate fits for a = +0.998, 0, and —0.998, to explore the dependence 
on the spin. A slightly better fit was obtained for the case of negative spin, with 
and ~5 higher with respect to the cases with zero or positive rotation. Thus, we conclude 
again that the quality of the current Suzaku data of 3G 382 does not allow us to constrain 
the sign of the black hole spin. Future deeper observations (>300 ks) are needed to this end. 
However, given that the fit is slightly better for a spin of —0.998 and that the overall model 
parameters are still consistent at the 90% whatever spin is assumed, in the following fits and 
in Table 4 we fix the spin to this value. The inclusion of the ionized refiection component 
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improves the fit by A^^ = 65 for four additional model parameters, which correspond to a 
detection probability of >99.99%. 

The estimated ionization of the disk reflecting material is log^ = 3.041;q;q4 erg cm, 
consistent with the Fe ion population being dominated by Fe XXV, with the Hea at E~6.7 keV 
being the most intense associated resonant hne (e.g., Ross, Fabian & Young 1999; Ross & 
Fabian 2005; Garcia & Kallman 2010). The fitted outer radius is rout='^Q-4 i"g and the in- 
clination i=29° ± 2°, consistent with the more simplistic fits using diskline or relline. The 
fit is very good, with a reduced chi-squared very close to unity, x^/dof = 4430/4168. The 
best-fit parameters of this broad-band fit are reported in Table 4, while the data, best-fit 
model and spectral ratios are shown in the left panels of Figure 4. 

5.5. Two ionized reflection components in 3C 382? 

It should be noted that the modeling of the soft excess with a simple black body com- 
ponent in the previous sections, besides providing a good representation of the data, might 
be difficult to explain from a physical point of view. In particular, from a systematic X-ray 
spectral analysis of PG quasars, Gierlinski & Done (2004) demonstrated that if the soft 
excess is modeled with a simple black body component, the resultant temperature is in the 
narrow range /cT~100-200 eV, despite the large range in Eddington luminosities and black 
hole masses of the sources. The temperature is also too high to be directly associated with 
standard Shakura-Sunyaev accretion disks. This suggests a possible atomic origin for the 
soft excess and a promising alternative physical explanation is in terms of a blurred ionized 
refiection component (e.g.. Crummy et al. 2006). However, the modeling of the soft excess 
is still required after the inclusion of the ionized disk refiection component in the previous 
section. When parameterized with a black body, it still provides a high improvement of 
the fit, Ax^=538 for two model parameters, which corresponds to a detection probabihty 
>99.99%. 

The presence of an additional mildly ionized refiection component is required to self con- 
sistently model the Fe Ka/K/3, Ni Ka emission lines and refiection hump. In particular, the 
energy of the Fe Ka is ~6.4 keV for a wide range of ionization, from Fe I to Fe XVII (Kallman 
et al. 2004). This gives rise to the possibihty of simultaneously modehng the soft excess, Fe K 
lines and reflection hump with an additional reflection component. Therefore, we performed 
a test replacing the previous pexrav, black body and Fe K Gaussian emission lines with a 
cutoff power-law continuum and an additional xillver reflection component. To account for 
the velocity broadening, we convolved the latter with a Gaussian profile with the same width 
of the Fe Ka emission line, i.e., cr=90 eV, which corresponds to FWHM~10,000 km/s (see 
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Table 4). 

The best-fit values are reported in Table 5. The power- law continuum, warm absorber 
and relativistic highly ionized reflection parameters are essentially unchanged. The ionization 
parameter of the mildly ionized reflection is log^=1.5 ± 0.03 erg cm. This is consistent 
with Fe XIII-XVII being the most abundant Fe ions (Kallman et al. 2004; Garcia, Kallman 
& Mushotzky 2011). As shown in the right panels of Figure 4, this provides a very good 
fit, with x^=4467 for 4174 degrees of freedom. The fit is only slightly worse than that with 
the model in Table 4. The reduced is 1-06 and 1.07 using the previous and present 
model, respectively. However, the present model has the strong advantage to physically 
self-consistently explain all the main spectral features in the broad-band spectrum with the 
lowest number of model components. This result is remarkably similar to what observed 
for local Seyfert galaxies (e.g., Nandra et al. 2007; de la Calle Perez et al. 2010; Patrick et 
al. 2010). 

The reflection fraction R is conventionally used to estimate the relative contribution of 
the reflection component with respect to the power-law continuum emission. This parameter 
is directly derived when using the pexrav model in XSPEC, as in Table 3 and Table 4. This is 
related to the solid angle covered by the reflecting gas as R—fl /27r. However, anisotropics, 
obscuration and relativistic effects close to the black hole can all influence this simple con- 
ventional relation for AGN, to the point that it can become physically meaningless (Crummy 
et al. 2006; Garcia, Kallman & Mushotzky 2011). Therefore, we define a phcnomenological 
reflection flux fraction, Rp, as the ratio between the unabsorbed 0.5-100 keV fiux of the 
refiection component and continuum. Its value is still comparable to the conventional pa- 
rameter, Rp^R (Garcia, Kallman & Mushotzky 2011). As reported in Table 5, this fraction 
is about 0.1 for both highly and mildly ionized reflection components, which means that each 
of them contributes ~10% on the overall broad-band spectrum. In particular, the mildly 
(highly) ionized component contributes ~20% (~30%) in the 0.5-2 keV, ~5% (~10%) in 
the 2-10 keV and ~10% (~5%) in the 10-200 keV with respect to the power-law continuum 
(see Figure 4 right, upper panel). 

Due to the wide-band coverage allowed by the BAT, very good constraints can be 
obtained on the power-law continuum energy cutoff. The latter value of Ec— 180 keV is 
completely consistent for both models presented in this section and § 5.4, see Table 4 and 
Table 5. This clearly suggests a thermal disk-corona origin for the primary power-law. Thus, 
it is unlikely that the X-ray power-law continuum of 3G 382 is due to a jet. 

In conclusion, the broad-band 0.6-200 keV spectrum of 3G 382 can be described by a 
very complex model composed by a Galactic absorbed power-law continuum, a black body 
soft excess, a warm absorber, an Fe Ka emission line and related Fe K/3, a Ni Ka emis- 
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sion line, a mildly ionized reflection component and a relativistic smeared ionized reflection 
component. Alternatively, it can be well represented also by a Galactic absorbed cutoff 
power-law continuum, a warm absorber and two ionized reflection components: a mildly 
ionized one which simultaneously reproduces the narrow Fe Ka emission line, soft excess 
and high energy hump and a highly ionized one which reproduces the broad relativistic Fe K 
line. These two models are indistinguishable from a statistical point of view, but the latter 
provides a more physically self-consistent representation of the data, with the lowest number 
of model components. 

6. The XMM-Newton Spectrum 

The spectral analysis of the EPIC pn observations was carried out using the heasoft 
V. 6.5.1 package. As shown in Figure 5, the overall 0.6-10 keV pn spectrum is very similar to 
the XIS one, albeit at lower S/N due to the lower exposure (factor 6) which compensates for 
the larger EPIC effective area compared to the XIS cameras. The ratios against a Galactic 
absorbed power-law continuum are reported in the middle panel of Figure 5. As we can 
see, the fit is very poor {x^/dof — 4149/1454) and the power-law photon index is P ~ 2. 
Intense emission residuals are present at low and high energes, similar to the Suzaku ones in 
the upper panel of Figure 2. These deviations indicate the presence of a soft excess and a 
reflection components. 

Because of the lower quality of the data and the spectral similarities, and because the 
source was in similar intensity states at both epochs (Table 2), we assumed for the EPIC 
pn the best-fit model previously derived for the Suzaku data from the fits to the broad-band 
spectrum discussed in § 5.2 and reported in Table 3. We assumed only a cold refiection 
component modeled with pexrav and fixed the cut-off energy and refiection fraction to their 
best-fit values. We also assumed only the presence of the narrow Fe Ka/K/3 emission fines, 
as indicated by the data. 

The fit to the EPIC pn data with this model is good, x^/rfo/=1567/1448. The best-fit 
model and spectral ratios are reported in the upper and lower panels of Figure 5, respec- 
tively. The power-law continuum photon index is r=1.87 ± 0.01. The EPIC data confirm 
the presence of a strong soft excess, with a blackbody temperature kT—lA2 ± 5 eV and 
normafization (1.72 ±0.03) x 10~^. The warm absorber is also detected, with column density 
A^i^=1.8to 2 X 10^^ cm~^ and ionization parameter log^=2.35 ± 0.05 erg s"-*^ cm, consistent 
with Suzaku (§ 5.2). 

An unresolved ((7=100 eV) Fe Ka line is detected at E=6.4 ± 0.03 keV, with a flux 
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7=(2.5 ± 0.4) X 10-^ ph s-^ cm-^ and EW=54l?^ eV. Due to the lower S/N of the pn data, 
there is no clear evidence for the associated Fe K/3 and Ni lines present in the XIS data, or for 
a broad, relativistic Fe K hne at 6.67 keV. The 90% upper limit for the latter is EW<130 eV. 

In summary, an archival 20 ks XMM-Newton EPIC pn observation of 3C 382 was ana- 
lyzed and found consistent with the best-fit model to the Suzaku XIS data. The observed and 
intrinsic (corrected from any intervening absorption) fluxes and luminosities in the 0.5-2 keV 
and 2-10 keV bands from the EPIC pn data are reported in Table 2. 

7. Summary of Observational Results 

This paper presented a 116 ks Suzaku observation of the Broad-Line Radio Galaxy 
3C 382 (^=0.057). A shorter (20 ks), unpubhshed XMM-Newton EPIC exposure was dis- 
cussed, as well as a Swift BAT spectrum from integrating 5-years of survey observations. 
The main results of the data analysis are: 

The soft X-ray spectrum - Below 2 keV, the XIS data indicate the presence of a warm ab- 
sorber, with parameters consistent with those derived from the analysis of a previous Chan- 
dra HETGS spectrum. The soft excess previously observed is confirmed by both Suzaku 
and XMM-Newton with an increase in flux between the two epochs of ~10%. This compo- 
nent can be well modeled by a phenomenological black body with A;T~100 keV, but a more 
physical parameterization suggests a possible origin from mildly ionized reflecting material. 

The Fe K Region - An intense narrow Fe Ka emission line is detected at E~6.4 keV, 
together with the associated Fe K/3 component. An additional weak unresolved emission 
line at E~7.5 keV is tentatively observed, possibly identifiable with Ni Ka. However, more 
interestingly, a prominent Fe K line with a relativistic profile is found; its center energy is 
consistent at 90% confidence with emission from ionized Fe, either or both FeXXV/FeXXVI 
in the energy range 6.55-7.10 keV. The profile is consistent with emission from a disk annulus 
between inner and outer radii rj„ ~ 10 and rout ~ 20 r^, respectively, for an assumed disk 
emissivity ^ — —2.5. An ionized disk refiection model was included to parameterize this 
component in the broad-band spectrum. 

The spectrum at energies >10 keV - The data at medium and hard X-rays are well described 
by a power law with F ~ 1.8, cutting off around Ec ~ 200 keV. Above 10 keV, contributions 
from both mildly and highly ionized reflection are present. Each of them contribute ~10% 
on the overall broad-band continuum, but the mildly ionized one has an higher impact at 
E>10 keV. 
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8. Discussion 

This paper presents new results from the analysis of our Suzaku as well as archival XMM- 
Newton and Swift BAT observations of 3C 382, one of the "classical", best-studied BLRGs 
from optical and radio samples. The combination of the Suzaku and BAT data provides 
an unprecedented sensitive coverage in the energy range 0.6-200 keV, useful to disentangle 
the various contributions to the total nuclear X-ray emission. At first sight, the broad-band 
X-ray spectrum of 3C 382 is remarkably Seyfert-like with most, if not all, the characteristics 
of a typical nearby radio-quiet AGN: a F ~ 1.8 power-law continuum, a warm absorber, soft 
excess, relativistic Fe Ka line, and a reflection bump above 10 keV. While typically radio- 
loud in other wavelengths, with giant radio lobes and a radio-to-optical Spectral Energy 
Distribution typical of other BLRGs (e.g., Grandi et al. 2001), 3C 382 appears exclusively 
radio-quiet from an X-ray spectroscopy perspective. 

The intrinsic photon index we measure for 3C 382 after accounting for the reflection 
components is F~1.8, similar to most radio-quiet Seyfert Is observed with Suzaku (e.g., 
Patrick et al. 2010). Previous observations of 3C 382 at X-rays (§ 2) established that the 
photon index varies in unison with the flux, in the sense of a softer spectrum for increasing 
X-ray intensity (Gliozzi et al. 2007 and references therein), a trend 3C 382 shares with 
radio-quiet, non-jetted AGN. Together with the constraints on the high-energy cutoff of 
the primary power law, this implies that thermal Comptonization dominates the emission 
below 100 keV in 3C 382, supporting the idea that the bulk of the X-ray continuum does 
not originate in a jet. At present, only an upper limit is reported to the gamma-ray GeV 
emission at the location of 3C 382 in the Fermi 15-months catalog (Abdo et al. 2010). This 
suggests that the contribution of a jet to the higher energies is weak, or occurs with a small 
duty cycle. Future Fermi monitoring will be necessary to shed light on the jet activity at 
high energies from the core of 3C 382. 

The high quality and low background of the XIS data allows us to model with high 
fidelity the Fe-K emission. We confirm the presence of a narrow line at an energy of ~6.4 keV, 
identified with the narrow component of the Fe Ka line. A similar narrow component, 
consistent with emission from the Broad Line Region (BLR), was detected in 3C 390.3 (S09) 
and 3C 120 (Kataoka et al. 2007), two of a handful of well-known, nearby BLRGs observed 
with Suzaku and is regularly observed in Seyfert Is (Patrick et al. 2010; Nandra et al. 
2007). This line can be modeled with a mildly ionized reflection component with log^~ 
1.5 erg s~^ cm (Table 5), which can also simultaneously explain the observed soft excess and 
part of the reflection hump. This contributes ~10% on the power-law continuum. 

Prom the width of the line of FWHM~10,000 km/s and assuming that the material is 
in Keplerian motion in the gravitational potential well of the super-massive black hole, with 
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estimated mass Mbh ^ 10^ M© (Marchesini et al. 2004), we can derive a typical distance of 
~0.3 pc. The similarity of the Fe Ka line width with that of the broad Ha line in the optical 
spectrum (Eraclcous & Halpern 1994) and the distance on sub-pc scales possibly suggest an 
origin of this reflecting material from the BLR. Using the distance estimate in the dcflnition 
of the ionization parameter ^ = Ljo„/nr^ (Tarter, Tucker & Salpeter 1969), where n is the 
average absorber number density and Lion is the source X-ray ionizing luminosity integrated 
between 1 Ryd and 1000 Ryd (1 Ryd=13.6 eV), we can derive the density of the material at 
that location. The absorption corrected ionizing luminosity of 3C 382 is Ljo„=8x 10^^ erg s~^. 
Substituting this value, we can estimate the density of the material of ~5 x 10^ cm~^. 

Moreover, we also detected a weak emission line at E~7.5 keV, most probably ascribable 
to Ni Ka fluorescence emission line from the same material of the Fe Ka. In particular, the 
Fe Ka is expected to be more intense than the Ni Ka due to the higher cosmic abundance 
of iron with respect to nickel, about ~20 times (Molendi, Bianchi & Matt 2003; Yaqoob 
& Murphy 2011). However, the ratios of the intensities of the two emission lines suggest a 
possible slight overabundance of nickel and/or an underabundance of iron of a factor of ~2 
with respect to the standard solar values (e.g., Anders & Grevesse 1989; Grevesse et al. 1996). 
Unfortunately, we note that the treatment of Ni lines is not included in the present version 
of the xillver code (March 2010). However, the possible slight underabundance of iron or 
overabundance of Ni has neghgible effects, within the la errors, on the parameters of the 
mildly ionized reflection component. 

Our most interesting result from the Fc K region modeling is the detection of a rcla- 
tivistic profile for the broad component of the Fe Ka line, confirming a previous tentative 
ASCA result (Reynolds 1997; S99). The Suzaku data quahty is sufficient to provide tight 
constraints on the disk parameters, albeit in a model-dependent way, in particular on the 
inner and outer radii of the emitting region and the disk inclination. We flnd that the bulk of 
the reprocessing in 3C 382 occurs in a small disk annulus, Ar ~ lOr^, located at a distance 
Tin ~ lOr^ from the central black hole and at an inclination i~25°-30°. Substituting the 
average distance of r~15rg in the ionization parameter equation and the estimated value 
of log.^^3 erg s~^ cm (Table 5), we can derive a density for the associated highly ionized 
reflection component of ~10^^ cm~^. This is a realistic estimate for the inner parts of an 
accretion disk in AGN (Garcia, Kallman & Mushotzky 2011). 

A similar finding - that the innermost disk regions do not contribute to the Fe K emission 
- was previously reported for other broad-lined radio-loud AGN, the BLRGs 3C 120 and 
3C 390.3 (Kataoka et al. 2007, S09), and the quasar/BLRG 4C+74.26 (Larsson et al. 2008). 
In these sources the inner radius of the Fe K emission was determined to be around lO-SOr^ 
from fits to the Suzaku data. Thus, a pattern is starting to emerge, whereas the innermost 
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regions of the accretion disk in broad-lined radio-loud AGN do not contribute to the Fe K 
emission. 

A number of plausible physical scenarios have been put forth to explain the above lack of 
Fe K reprocessing in the innermost radii of the BLRG disks (see S09 and discussion therein) . 
These include a highly ionized ion torus/ ADAF occupying the inner disk, obscuration by 

the base of a jet, or lack of illumination of the inner radii of the disk by the primary power 
law, due to a mildly beamed source of continuum. An additional possibility is offered in the 
context of the "flux-trap" scenario (Reynolds et al. 2006; Garofalo 2009). 

It is generally believed that jet power is, at least in part, linked to the spin of the central 
black hole. Radio- loud AGNs have been claimed to harbor rapidly spinning black holes (e.g., 
Nemmen et al. 2007, and references therein). A scenario for jet power production is provided 
by the "flux trap" model (Reynolds et al. 2006; Garofalo 2009), which suggests that the 
plunging region between accretion disks and black holes is fundamental in producing strong, 
spin dependent, horizon-threading magnetic fields, and therefore powerful jets. Noteworthy, 
it has been demonstrated that the Blandford & Znajek (1977; BZ) mechanisms is maximized 
for black hole spin a ~ — 1 (retrograde) and stronger jets/outfiows are expected in this case 
(Garofalo 2009). This is a direct consequence of the fact that the radius of marginal stabihty 
for retrograde black holes, rjsco 9rc, for a = —0.998, is much larger than in prograde, 
f^isco ~ 1.23rg for a = +0.998, and therefore the plunging region is more extendend and 
the fliix-trapping effect is magnified. In principle, this could be directly tested with X-ray 
spectroscopy because in the two cases the relativistic Fe K line profile changes significantly 
(Garofalo, Evans, & Sambruna 2010). 

As discussed in § 5.3 and § 5.4, we tested the retrograde/prograde scenario by performing 
spectral fits to the data with a model that allows the black hole spin to assume negative 
values (Dauser et al. 2010). While the fit is statistically satisfactory, unfortunately the spin 
parameter is unconstrained due to the insufficient S/N ratio of the data. However, because 
of the shape of its Fe K line, 3C 382 quahfies as the best candidate so far for determining 
the black hole spin in a radio-loud AGN through Fe K line spectroscopy. Longer exposures 
with Suzaku are required. 

3C 382 shares with two other BLRGs (3C 390.3 and 3C 120) the evidence in its X- 
ray spectrum for emission from an ionized disk, with similar ionization parameters ^ ~ 
2300 — 2700. This supports the previous suggestion by Ballantyne et al. (2002) that ionized 
disks play a significant role for the X-ray emission of powerful radio-loud AGN. Evidence for 

ionized material in the nuclear regions of radio-loud AGN is also provided by the detection 
of extremely ionized outfiows/winds with velocities >10,000 km/s, and properties similar to 
radio-quiet AGN outfiows (Tombesi et al. 2010a,b; see also below). 
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Finally, we compare the X-ray properties of 3C 382 and other BLRGs observed with 
Suzaku to Seyfert Is. We started this comparison in our previous paper on 3C 390.3 (S09), 
where we outhned the average BLRG properties - X-ray photon index, Fe K hne widths, and 
cold reflection strengths - with the then available most current X-ray data for Seyferts from 
XMM-Newton EPIC observations (Nandra et al. 2007). We noted how the XMM-Newton 
data for radio-quiet AGN suggested a much larger spread of reflection fractions and line 
widths for Seyferts than in the ASCA and RXTE era, and how this spread lead to a blur in 
the distinction between radio-loud and radio-quiet sources. In particular, in BLRGs X-ray 
reflection and continuum properties occupy one end of the distribution for Seyferts, with 
significant overlap. We concluded that the division between radio-loud and radio-quiet AGN 
was blurred, as far as their X-ray spectral properties were concerned (S09). 

A more proper comparison is now possible thanks to the recent publication of Suzaku 
and Swift BAT observations for a relatively large (6 sources) sample of Seyfert Is (Patrick et 
al. 2010). From an accurate modeling of the 0.6-100 keV spectra using both cold and ionized 
reflection models, Patrick et al. (2010) found broad relativistic profiles for the Fe K emission 
lines in most Seyferts, with average equivalent width for the sample of EW=119±19 eV and 
ionized emission lines at 6.7 and 6.97 keV being relatively common. Intriguingly, the Fe K 
emission seems to arise from tens of r^ from the central black hole (see also Nandra et al. 
2007); fits with models allowing the black hole spin to vary, although in the 0-1 range only, 
yield typically intermediate values, a ~ 0.7. Both cold and ionized refiection is required 
to describe the data above 10 keV; the cold refiection fraction, however, varies from source 
to source from less than 0.4 up to 2.5 (see Table 4 in Patrick et al. 2010). The intrinsic 
X-ray photon index is Tx ~ 1-8 — 2.2, consistent with Comptonization models. However, 
the ionization parameter of the ionized disk reflection seems to be lower in the Seyferts, 
consistent with neutral or lowly ionized Fe, while interestingly so far the BLRGs seem to 
suggest high ionization disk lines (He and H-like iron), although based on small number 
statistics so far. 

Thus, comparison of the broad-band 0.6-200 keV spectra of BLRGs and Seyfert Is 
obtained with Suzaku and the BAT conflrm and reinforce the conclusion (S09) that there is 
a continuum of X-ray properties between the two AGN subclasses, with radio-loud sources 
clustered at one end of the distribution. An important consideration, however, is that in 
both subclasses - aside for a few egregious examples - the innermost ~ lOr^ of the accretion 
disk do not contribute signiflcantly to Fe K emission. The reason for this may be different in 
the two subclasses, and would remain of central importance for the radio-loud/radio-quiet 
AGN division. 

Another major progress since S09 is the discovery that radio-loud AGNs have ionized. 
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outflowing material from the central nuclei, similar to Seyferts (e.g., Tombesi et al. 2010a,b), 
contrary to what stated in S09 based on the then available evidence. Indeed, our Suzaku ob- 
servations of 5 BLRGs provided evidence for ultra-fast disk outflows in 3/5 sources (Tombesi 
et al. 2010b) with velocities Vout ~ 0.1c, kinetic power ~ \{)'^^-^'^ g-i^ ^nd mass outflow 
rates comparable to the accretion rates. As mentioned above and in § 2, 3C 382 itself presents 
evidence for a highly ionized, kpc-scale outflow through a series of soft X-ray absorption fea- 
tures detected with the Chandra and XMM-Newton gratings (Reeves et al. 2009; Torresi et 
al. 2010), although with a much lower velocity of ~1,000 km/s. The presence of ionized 
material along the line of sight is puzzling and its relation with the relativistic jet must be 
studied in much detail. As already discussed in Tombesi et al. (2010b), we did not detect 
any significant highly ionized Fe K absorption line in the present observation of 3C 382. 

The overlapping X-ray continuum and refiection properties of BLRGs and Seyferts point 
to a common accretion fiow structure, with similar radiative processes in place that account 
broadly for the X-ray emission. The reason why powerful, relativistic jets are present in 
BLRGs but not in Seyferts with similar engines could very well be identified in a process 
unrelated to the accretion dynamics - i.e., the black hole spin. Several theoretical studies 
along this line are underway, and the wealth of high-quality data that Suzaku continues to 
provide is essential to test the model predictions in the X-ray band, paving the way for the 
upcoming high-precision spectroscopy with the Astro-H observatory. 

This research has made use of data obtained from the High Energy Astrophysics Sci- 
ence Archive Research Center (HEASARC), provided by NASA's Goddard Space Flight 
Center, and of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, under contract with the National 
Aeronautics and Space Administration. RMS and FT acknowledge support from NASA 
through the Suzaku and ADAP/LTSA programs. FT thank T. Dauser for help in using the 
relline model. FT thank J. Garcia for help in using the xillver model. LB acknowledges 
financial support from the Spanish Ministry of Science and Innovation through a "Juan de 
la Cierva" fellowship, research grant AYA2009-08059. The authors thank the referee for 
suggestions that led to important improvements in the paper. 



A. Background and Calibration tests for the XIS 

We performed detailed background and calibration tests for the three XIS detectors 
to exclude possible issues affecting the Fe K region, and thus the results of our modeling 
described above. Background contamination could alter the profile of the Fe line; however. 
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for the 3C 382 observation analyzed here, we estimate that the background contribution to 
0.6-10 keV is neghgible, amounting to <4-8% of the source counts. 

It is well known that the XIS CCDs contain a ^^Fe calibration source, located on two cor- 
ners of the XIS chips, which produces a characteristic X-ray hne from Mn Ka at E=5.895 keV. 
Its intrinsic width is a"<50 eV and it is used to calibrate the energy and width of the iron 
line (Koyama ct al. 2007). Wc can clearly exclude any possible significant contamination of 
the broad line detected at the rest-frame energy of E~5.8 keV from this calibration source 
on board the XIS. In fact, we checked the images of the three separate XIS cameras and find 
that the source is located in the center of the CCD, far from the calibration sources, which 
are in two corners. Moreover, the observed energy of the broad line is at E=5.45 keV, which 
is not consistent with that expected from the calibration source, and the calibration line is 
intrinsically narrower compared with the ~400 eV of the 3C 382 feature. 

As already discussed in § 3.1 the combined XIS-FI and XIS-BI fit is supported by the 
consistency of the power-law slope and 2-10 keV flux within ~2%. Moreover, as already 
reported in § 4, we checked that the presence of two low amplitude (~5%) spikes in the 
XIS 3 light curve, see Figure 1, have neghgible effects on the overall spectral analysis. This 
is expected from the fact that we focused on the average spectrum and the weak short time- 
scale spikes involved less than 20% of the exposure. These are most probably ascribable 
to weak instrumental fluctuations. Wc directly analyzed the three XIS cameras separately 
and ensured that there is an agreement at ~1% for both the power-law continuum slope 
and 2-10 kcV flux between the XIS and XIS 3 and at ~2% between the front illuminated 
XIS 0-3 and back illuminated XIS 1. 

Finally, the spectra of the individual XIS camera were fitted separately for a consistency 
check of the emission line parameters. Unfortunately, the quality of the individual camera 
spectra does not allow us to leave all line parameters free to vary during the fit. Using the 
best-fit model determined from the analysis of the joint spectra (see Table 3), and freezing 
the fines energies and widths, we find that the EWs of the emission fines are always consistent 
at 90% among the separate XIS instruments and the combined XIS-FI and XIS-BI fit. 
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Table 1. 3C 382 observations log. 



Satellite 


OBSID 


Date 


Instrument 


Net expo 
(ks) 


Rate 
(cts/s) 


Suzaku 


702125010 


2007 April 27 


XIS-FI 


116 


2.399/1.219 








XIS-BI 


116 


1.613/1.177 








PIN 


113 


0.114 


XMM-Newton 


0506120101 


2008 April 28 


EPIC pn 


21 


15.650/4.220 


Swift 


J1835.0+3240 




BAT 


4276 


0.001 



Note. — The count rates refer to each instrument separately and are in the 0.5-2/2-10 keV 
band for XIS-FI, XIS-BI and EPIC pn, 14-70 keV band for the PIN and 15-100 keV band 
for the BAT. 



Table 2. Flux and luminosities from Suzaku and XMM. 



Energy 


Suzaku 


XMM 


(keV) 


Flux Lum 


Flux Lum 


0.5-2 


2.42(3.10) ±0.02 1.92(2.51) ±0.02 


2.59(3.44) ± 0.02 2.05(2.79) ± 0.01 


2-10 


4.09(4.12) ±0.01 3.25(3.28) ± 0.01 


3.83(3.88) ±0.02 3.05(3.09) ± 0.01 


15-200 


8.10(8.10) ±0.40 6.50(6.53) ± 0.30 





Note. — The observed (un-absorbed) flux is in units of 10 erg s ^ cm ^. The 
observed (un-absorbed) luminosity is in units of 10^** erg/s. 
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Table 3. Best-fit of the E=0. 6-200 keV broad-band spectrum. 



Pexrav 



r 


Ec 
(keV) 


R 


i 

(deg) 




1.859 ± 0.004 




0.36 ± 0.06 


30 


4416/4167 


Warm absorber 


(1020 cm-2) 


log? 

(erg s^-' cm~2) 


Vout 

(km/s) 






5.8 ± 1.5 


^■^^-0.04 


1000 






Black body 


feT 
(eV) 


Norm 
(10-4)^ 








100 ±4 


1.4 ±0.2 








Gaussians 


E 
(keV) 


cr 
(eV) 


/ 

(10-5 ph s-l cm-2) 


EW 
(eV) 


ID 


6.41(6.06) ±0.01 

= 7.06(6.67) 
7.51(7.10) ± 0.03 
5.70(5.39) ± 0.09 
6.91(6.53) ±0.02 


= 10 
= 10 


3.3 ±0.2 

0.7 ±0.1 
3.0 ±0.4 
0.7 ±0.1 


68 ±9 

10 
18 ± 7 
55± 12 
17 ±6 


Fc Ka 
Fe K/3 
Ni Ka 



Note. — Energy of the line, rest-frame (observed frame). Errors are at the Icr level. 

^Normalization of the blackbody model in units of Lag/Djg, where Lag is the source 
luminosity in 10^^ erg s"^ and Dio the distance of the source in units of 10 kpc. 
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Table 4. Best-fit of the E=0. 6-200 keV broad-band spectrum, including the ionized 

reflection component. 



Pexrav 


r 


(kcV) 


R 


i 

(dcg) 




i.Vil = 0. ()():-! 


1S7 1 :!1 


O.ir, ± ().()-, 


:-!() 


1 !:-!()/ il()8 


Warm absorber 


(1020 cm-2) 


log? 
(erg cm~^) 


Vout 

(km/s) 






8.5 ±1.5 


^■"'-'-0.06 


1000 






Black body 


fcT 

(cV) 


Norm 

(10-*)=^ 








104 ±3 


1.4 ±0.5 








Gaussians 


E 
(keV) 


(eV) 


/ 

(10-5 ph s-1 cm-2) 


EW 
(eV) 


ID 


6.41(6.06) ±0.01 

= 7.06(6.67) 
7.51(7.1) ±0.03 


= 10 


2.1 ±0.2 
0.5 ±0.1 


42 ±6 
6 

13 ±7 


Fe Ka 
Fe K/3 
Ni Ka 


Relconv 


13 


a 


fout 


i 

(deg) 




-2.5 


-0.998 


26t5 


29 ±2 




Xilh-er 


log? 
(erg cm) 


Norm 
(10-8) 










^■'-0.3 









Note. — Energy of the line, rest- frame (observed frame) . Errors are at the 1<t level. 

''Normalization of the blackbody model in units of L39 /D\q, where Lag is the source 
luminosity in 10^^ erg and Dio the distance of the source in units of 10 kpc. 
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Table 5. Best-fit of the E=0. 6-200 keV broad-band spectrum, including two ionized 

reflection components. 



Cutoff power-law 


r 


Ec 
(kcV) 






1.737 ± 0.003 


175+^* 




4467/4174 


Wcirm absorber 


(1020 cm-2) 


log? 
(erg cm~2) 


Vout 

(km/s) 




8.3 ±1.5 


''•^^-o.oe 


1000 




Xillver (mildly ionized) 


log? 
(erg cm) 


Norm 
(10-6) 


Rf 




1.54 ±0.03 


3.0 ±0.3 


0.10 




l-i,t;icoii\' 


/3 


a 


'^out 


i 

(deg) 


-2.5 


-0.998 


25l« 


30 ±3 


Xillver (highly ionized) 


log? 

(erg s^-"- cm) 


Norm 
(10-8) 


Rf 




2.93 ± 0.04 


4-^-0.3 


0.11 





Note. — Energy of the line, rest-frame (observed frame). Er- 
rors are at the Icr level. 

^Normalization of the blackbody model in units of L:it)/Dff^, 
where L39 is the source luminosity in 10^^ erg s-^ and Dio the 
distance of the source in units of 10 kpc. 
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Fig. 1. — Light curves of 3C 382 at soft (a, left) and hard (b, right) X-rays from the Suzaku 
XISO-3 observations. The data from the three detectors are plotted in separate panels for 
clarity. Bins are 5760 sec, or 1 orbit. 
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Fig. 2.— Combined XIS-FI (black), XIS-BI (red), PIN (green) and BAT (blue) broad-band 
spectra in the E=0. 6-200 keV energy band. Data have been rebinned to S/N=40 for XIS and 
S/N=12 for PIN only for plotting. Upper panel: ratios with respect to a Galactic absorbed 
power-law continuum. The Fe K band between 5.5-7.5 keV is excluded. Middle panel: ratios 
after the inclusion of a blackbody soft excess and a cold reflection component. Lower panel: 
ratios after the inclusion of a warm absorber and five Gaussian emission lines in the Fe K 
band. 
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Fig. 3.— Zoomed spectrum of 3C 382 in the Fe K band, E=4-9 keV. The XIS-FI (black 
sohd hne) and XIS-BI (red dotted hne) have been both rebinned to 100 eV energy bins only 
for plotting. Upper panel: data/model ratios with respect to the best-fit continuum model 
without the inclusion of emission lines in the Fe K band. The arrows indicate the detected 
Gaussian emission lines. Middle panel: best-fit model composed by three narrow Gausssian 
emission lines and a diskline profile. Lower panel: data/model ratios with respect to the 
model including three Gaussian emission lines and a diskline. 
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Fig. 4.— Combined XIS-FI (black), XIS-BI (red), PIN (green) and BAT (blue) broad-band 
spectra in the E=0. 6-200 keV energy band. Data have been rebinned to S/N=40 for XIS 
and S/N=12 for PIN only for plotting. Left: broad-band model including the black body soft 
excess, Gaussian emission lines and cold pexrav reflection component, see § 5.4 and Table 4. 
Right: broad-band model including only two ionized reflection components, see § 5.5 and 
Table 5. Upper panels: data with best-fit model. Lower panels: ratio with respect to the 
best-fit model. 
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Fig. 5.— XMM-Newton EPIC pn spectrum of 3C 382 in the E=0.6-10 keV band. Data 
have been rebinned to a S/N of 20 only for plotting. Upper panel: data and best-fit model. 
Middle panel: ratio against a Galactic absorbed power-law continuum. Lower panel: ratio 
against the best-fit model. 



